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Abstract 

Ever since the first observations that we are Hving in an accelerating 
universe, it has been asked what dark energy is. There are various ex- 
planations all of which with have various draw backs or inconsistencies. 
Here we show that using a dark spinor field it is possible to have an equa- 
tion of state that crosses the phantom divide, becoming a dark phantom 
spinor which evolves into dark energy. This type of equation of state has 
been mildly favored by experimental data, however, in the past there were 
hardly any candidate theories that satisfied this crossing without creating 
ghosts or causing a singularity which results in the universe essentially 
ripping. The dark spinor model converges to dark energy in a reasonable 
time frame avoiding the big rip and without attaining negative kinetic 
energy as it crosses the phantom divide. 

1 Introduction 

An increasing number of independent observations indicates that we are living 
in an expanding universe where the expansion itself is accelerating [U [51 12] . 
It has been accepted that this requires some additional negative-pressure mat- 
ter source, named dark energy. The simplest model explaining this accelerated 
expansion is the cosmological constant A which corresponds to an unusual equa- 
tion of state w = P/ p = —1. Although the standard model of cosmology, the A 
cold dark matter (ACDM) model, fits the present data very well, the numerical 
value of the cosmological constant is about 120 orders of magnitude smaller 
than the vacuum expectation value predicted by quantum field theory. This 
smallness problem can be addressed by considering dynamical models, where 
the field slowly rolls down some potential where the effective equation of state 
Weff converges to w^s = ^1- In the past it has been believed that this value 
should be approached from above. However, there has been a recent interest in 
phantom models where the dark energy equation of state is approached from 
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below'u;< -1, see [1 [3 [3 [HI [IHl M [11 [II [II [111 [M [13 [IB IIS] • These 
models, as counter intuitive as they may appear, are not excluded by current 
data dig. 

Figure [J taken from [5] , shows data taken from the cluster abundance, 
supernovae, quasar-lensing statistics and the first acoustic peak in the cosmic 
microwave background radiation (CMB) power spectrum, which all imply that 
there is convergence to dark energy. However, when the parameter space is 
expanded to include the phantom region the data does not rule out a convergence 
from below [S]. In fact, the data even seems to indicate that this is quite likely. 

Considering such phantom energy models, the outcome of such a universe 
dominated by phantom energy is very different to any outcome we are accus- 
tomed to. The evolution parameter increases quicker than the horizon and 
quickly starts to pull apart gravitationally bound objects and then finally any 
bound objects, including the strong force! However, since cosmologists turned 
their attention from a constant equation of state for dark energy to one where it 
evolves through time it is believed that with or without crossing the divide the 
equation of state should eventually converge to w = — 1 due to the success of 
the ACDM model, and that observationally it is so close to this value at present. 
The majority of such dynamical dark energy models is based on evolving scalar 
fields with a suitably chosen potential. However, a canonical scalar field can- 
not cross the phantom divide and therefore such models have to modified, for 
instance by considering fields with negative kinetic energy. 

Within theories of modified gravity one is left with the option of either chang- 
ing the geometrical side or the matter side of the gravitational field equations, 
changing either gravity itself or altering the matter content which contributes 
to the forces, respectively. Those two approaches are not entirely independent 
as many modified theories bring the new geometrical quantities to the matter 
side, ultimately changing the energy content of the universe, however allowing 
for a new interpretation. 

Here, we follow a different approach, namely we model the evolving dark 
energy component of the universe to be of the form of a single dark spinor field. 
Dark spinors, originally dubbed the Elko spinor [5D], are similar to Majorana 
spinors but acquire the full four degrees of freedom of a Dirac spinor due to 
their helicity structure. They couple to the Higgs mechanism via XXH^H and 
weakly to the electromagnetic field via Xl^y"' ,j^]XFab, however in the latter case 
this coupling is heavily constrained because of the masslessness of the photon, 
making them a candidate for dark matter. The idea of one field explaining both 
dark matter and dark energy has already been discussed in various approaches, 
see e.g. [2T1[T8]. One possible mass range for our dark spinors is in the m ~ MeV 
range. 

The dark spinor field is based on the eigenspinors of the charge conjuga- 
tion operator, for a general introduction we refer the reader to [5D]. For the 
subsequent developments of dark spinors as an alternative model to scalar field 
inflation, see [22l [23l [24]. For instance, in [25] it could be shown that dark 
spinors have the effect of suppressing the low multipole amplitude of the pri- 
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Figure 1: Current constraints to the w — flm parameter space. The red solid 
curves show the age (in Gyr) of the Universe today (assuming a Hubble parame- 
ter Ho = 70kmsec~^ Mpc~^). The light shaded regions are those allowed (at 2a 
confidence level) by the observed cluster abundance and by current supernova 
measurements of the expansion history. The dark orange shaded region shows 
the intersection of the cluster-abundance and supernova curves, additionally re- 
stricted (at 2f7 confidence level) by the location of the first acoustic peak in the 
cosmic-microwave-background power spectrum and quasar-lensing statistics. 



mordial power spectrum, while at the same time to provide the usual inflationary 
features and predicting a small scale power spectrum in agreement with scalar 
field inflation, see in particular [26l [27l [28] for dark spinors in cosmology. 



2 Cosmological dark spinor field equations 

The standard model of cosmology is based upon the fiat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric 

ds^ =dt^ ~a{tf{dx^ +dy^ + dz^), (1) 

where a{t) is the scale factor and t is cosmological time. The dynamical behavior 
of the universe is determined by the cosmological field equations of general 
relativity 

Ra0 - -^Rgafl = JpTa0, (2) 

where Mp\ is the Planck mass which we use as the coupling constant, l/M^^ = 
SttG and c — I. Ta/3 denotes the stress-energy tensor, which for a homogeneous 
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and isotropic cosmology takes the form 

To,p = diag(/9, a^P, a^P, a^P). (3) 
The cosmological field equations are given by 

p + 3ff(p + P) = 0. (5) 

where the dot denotes differentiation with respect to time t and the Hubble 
parameter H is defined hy H = d/a. 

Let us consider a homogeneous single dark spinor field. Following [551 127] , 
the effective Lagrangian density of this field can be written in terms of the scalar 
field if as 

^ = + ^i/V^ - V{^). (6) 

If the potential V{(p) contains a standard mass term V{ip) — to^(^^/2, then we 
can rewrite the Lagrangian as 

^ = i^2 + ^i/V-^™V. (7) 

This allows us to interpret the explicit presence of the Hubble parameter in the 
action as an effective mass term where the mass changes as the universe evolves, 
and we have 



TO 



H'. (8) 



If the universe undergoes a phase of accelerated expansion, the Hubble param- 
eter is approximately constant. Depending on the ratio m/H, it is possible 
for models to attain a negative value for m'^g without creating ghost of having 
negative kinetic energy. This is possible because of the additional coupling of 
geometry to the matter field we are considering, a spinor opposed to a scalar 
field. 

The energy density and the pressure of the dark spinor field are given by 

= + V{^) - ^i/V, (9) 

= - V{^) + \h\\ (10) 

These two equations have the important property of leaving the acceleration 
equation unchanged 
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The spinor field's potential energy may yield an accelerated expansion of the 
universe. It should be noted that the energy density and the pressure now 
explicitely depend on the Hubble parameter. These additional terms are present 
because the covariant derivative when acting on the spinor field also acts on its 
spinorial part. One way of interpreting this 'coupling' in Eq. ^ is to regard 
the effective mass of the particle to depend on the Hubble parameter [33] and 
therefore on the evolution of the universe. Another interesting way of reading 
this equation is to regard the gravitational coupling to be time dependent |27] . 
The effective equation of state of the dark spinor field is given by 

and its more complicated form, when compared with the scalar field, also shows 
that a crossing of the phantom divide is possible without negative kinetic energy. 
If we for the moment take the canonical mass term for the potential and assume 
that </)^ <C y{}P), the effective equation of state becomes 

w^cff ^ ^-^Wi = - 1 \Ji 2 ^ 13 

TO^ — 1 — jH^ /m^ 

Therefore, we expect dark spinor models with canonical mass term to cross the 
phantom divide, however, we do not expect the effective equation of state to 
diverge since once the field is slow-rolling in its potential, there exists a lower 
bound to WcS- 

While there is a large number of potentials considered in the literature in 
the context of dynamical dark energy models, we will restrict our attention to 
power counting renormalizable potentials. As the dark spinor field has mass 
dimension one, the two allowed potentials are 

V^{^) ^ \ir?^\ (14) 

and 

V^{^) = \ir?p^ + )a^\ (15) 

where the V\ {(p) is the aforementioned canonical mass term and V2 {f) includes 
the self interaction term with a being a dimensionless coupling constant. 

3 Dark spinors as dark energy 

Let us start by solving Eq. ^ for the Hubble parameter, using Eq. ^ we find 



x/3Mpi ^1 + (^/Afpi)V8 
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Figure 2: Equation of state for Vi{(p): Mp\ — 1, (p(0) = 1 and w{0) = —1/3. 
With {0.002,0.001} = {red (higher), blue (lower)} 



Hence, using the first equation, the energy density of the dark spinor field is 

It is precisely this latter form of the energy density which motivated [57] to in- 
terpret (fTS]) as inducing a time-dependent gravitational coupling by considering 
Gtt = St^GcsP^p where p^p is the standard energy density of a scalar field. 



3.1 Phantom dark energy models 

Now, we consider the conservation equation ([5]) with (|16p and (jlSp and numeri- 
cally solve the resulting equation for Lp{t) and substitute into Eqs. ^ and (|10p 
to obtain the evolution of the effective equation of state Wcs — P/ P and plot it 
as a function of the evolution parameter a{t). We chose our initial conditions to 
be w[Q) — {1/3,0,-1/3,-2/3}, the first two representing radiation and dust, 
respectively. The initial conditions with w{Q) < —1/3 correspond to an initially 
accelerating universe. Small changes in these initial conditions do not alter the 
late-time asymptotic behavior of the solutions. 

Fig. [2] shows the dynamical behavior of the effective equation of state con- 
sidering the potential Vi with Afpi = 1, (p{0) = 1 and w{0) = —1/3. For the 
two different mass values it is possible to see that the effective equation of state 
almost immediately drops below the phantom divide. During the subsequent 
evolution, w begins to increase as further shown by Fig. [3] to the desired dark 
energy value. From Fig. [3] it is also evident that our model is practically indis- 
tinguishable from dark energy modeled by a cosmological constant, long before 
recombination when the scale factor is about a{t) — 10^"^. 
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Figure 3: Equation of state for V\{ip): — 1, (p(0) ~ 1 and ti;(0) = —1/3. 
With {0.002,0.001} = {red (higher), blue (lower)}. 
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Figure 4: Equation of state for Vi ((/?): A/pi = 1, (^(0) = 1 and u;(0) = 1/3. With 
= {0.002,0.001} = {red (higher), blue (lower)}. 



We obtained very similar results for other initial values of the equation of 
state: w(^) = 1/3, wifi) = and w{0) — —2/3. We have shown, for comparison, 
results from w(0) = 1/3 in Fig. [4] and Fig. [5l they qualitatively agree with the 
results presented in Figs. [2] and [3l respectively. 

3.2 Diverging models 

The next natural step is to consider a potential which allows for a self-interaction 
terms. As discussed above, the inclusion of this term also yields a power counting 
renormalizable field theory. 

When we include the self interaction term in the potential, ^2(95), we find, 
for all initial values of w and a — I, that the effective equation of state always 
diverges to —00, see Fig. ©. We have also checked, although not included 
in here, that our numerical results for V2{(p) converge to results for Vi{ip) as 
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Figure 5: Equation of state for Vi{(p): Mpi = 1, ip{0) = 1 and w{0) = 1/3. With 
= {0.002,0.001} = {red (higher), blue (lower)}. 
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Figure 6: Equation of state for V2(</?): Mpi — 1, w{0) = 1/3 and a = 1. With 
{4,0.02} = {red (higher), blue (lower)}. 



a 0. Although it might be possible to construct models with finely tuned 
initial conditions such that the divergence of the equation of state would happen 
in the future, such models are very unlikely. Hence, we are lead to conclude that 
a dynamical dark energy model based on dark spinors requires their potential to 
be of the simplest form, namely a canonical mass term, without self interaction. 

3.3 Oscillating models 

Lastly, we found another set of interesting results where the equation of state 
oscillated between w = 1 and w = — 1 for all time. The oscillation of the 
equation of state is very rapid, as can be seen in Fig. [3 making them unphysical 
as a dark energy model. This qualitative behavior does not change if we include 
the self-interaction term. However, if such a model could be modified it would 
be a prime candidate for models where the field changes its characteristic from 
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being dark matter at early times to become dark energy at late times, see 
also [291 [301 [311 [31 [33]. 
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Figure 7: Equation of state for potential ViiLp): = 1, to = 0.1 with 

initial conditions chosen such that u'(O) = {1/3,0,-1/3,-2/3}, respectively 
{blue (longdashed, red (medium dashed) , green (dashed), cyan (short dashed)} 



4 Conclusions and discussions 

A dark spinor field is able to provide a possible model for dark matter as it 
couples mainly via the Higgs mechanism, but has heavily constrained interac- 
tions with the electromagnetic field. Dark spinors have a predicted MeV mass 
range and therefore experimental predictions can be formulated and possibly 
measured at the LHC. Our results now show that the dark spinor field is also 
capable of having a dynamical equation of state which crosses the phantom di- 
vide and asymptotes to w = —1. This makes it a viable candidate for dark 
energy which cannot be ruled out experimentally. 

Unlike previous phantom models, dark spinors do not obtain negative kinetic 
energy on crossing the phantom divide, due to both p and P depending on the 
Hubble parameter and therefore these models do not create ghosts. According 
to |5] the equation of state must not stay below the divide but converge to dark 
energy, therefore the dark spinors' potential is of the simplest form, a canonical 
mass term to^i^^/2. Our dark spinor model does not require a modification 
of general relativity, leaving one of the most successful models in theoretical 
physics untouched. 

Due to the interesting nature of dark spinors, they have been shown to give 
other unique properties not found with other matter sources considered in the 
past. For now, in a cosmological setting, dark spinors are providing intriguing 
results in having the potential to be the best candidate dynamical dark energy 
model at hand. 
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